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Abstract 

A straightforward derivation of the effective quark Lagrangian is 
presented for the topologically neutral chiral broken phase of the dilute 
instanton gas. The resulting quark Lagrangian is a nonlocal NJL 
- type and contains 4q, 6q,... vertices for any number of flavours. 
Correspondence of this result with previously known in literature is 
discussed in detail. 

1 Introduction 

The instanton gas model has a long history [1,2,3]. Recently lattice measure- 
ments of quark and gluon correlators in coordinate space [4] have revealed 
important contributions of instantons, which gives an additional impetus to 
the model [5]. 

On the theoretical side the instantons have been proved to solve the U(l) 
problem [6], and to generate a specific effective quark Lagrangian (EQL) on 
a net topological charge [7]. 

In the topologically neutral system, such as the QCD vacuum, instantons 
have been proposed as a driving mechanism for the chiral symmetry breaking 
(CSB) [1,2,8]. The quantitative theory of light quarks (TLQ) in the instanton 
gas has been suggested in [9] and later elaborated in [10,11,12]. 
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There are a considerable number of papers which use the method for 
quantitative predictions and comparison with experiment (as an example see 
[13])- 

The essential ingredient of TLQ is the assumption of the dominant role of 
the quark zero modes on each instanton, and the resulting EQL is proposed 
in the form of the 2Nf x 2Nf determinant of quark fields with coefficients 
(vertices) M(p,p') expressed through quark zero modes. 

As a result TLQ predicts CSB and the creation of chiral quark mass at 
all instanton densities; after bosonization the effective bosonic Lagrangian 
was derived [9, 10] containing expected Goldstone bosons. 

In view of numerous successful applications we believe it is important to 
test the accuracy of the method. 

In its modern form [12] the TLQ contains some approximations, first of 
all the zero-mode approximations (ZMA) introduced already in [9], and the 
basic role in TLQ is played by the ansatz for EQL [10-12] of the 'tHooft form 
[7]. The ZMA was questioned in [14], and it was argued there that nonzero 
modes play an important role in addition to zero modes. 

In the present paper a new method is proposed, where all modes are 
accounted for automatically and we are avoiding standard approximations of 
TLQ. 

In the next section we demonstrate a straightforward and exact way to 
derive the EQL in the topologically neutral dilute instanton gas. Our results 
differ from the existing in several points and lay ground for a new formalism 
of light quarks in the instanton gas. 

First of all, our EQL contains a set of quark terms starting from 4q for 
any number of flavours. The 4q term differs from that of TLQ and is similar 
instead to the extended NJL type Lagrangian [15]. The structure of our EQL 
makes it possible to create chiral mass and CSB, and bosonization yields the 
Nambu-Goldstone boson Lagrangian. 

All coefficients in our EQL are explicitly known and this enables one to 
do quantitative predictions. In particular the appearance of the 4q vertex 
for any Nf is important from the point of view of theory (such as effective 
top-condensate model of electroweak interactions [16]) and applications (OZI 
violation vertex, quark-pair creation, strangeness mixture). 

In our formalism there are several ways to find the chiral quark mass 
entering the averaged quark propagator. One way exploits Dyson-Schwinger 
equations for our EQL in connection with the mean-field (Hartree) approx- 
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imation. This is equivalent to the direct summing of leading trms in the 
1/N C expansion. In the one-flavour case one obtains an integral equation for 
the chiral mass, which coincides with the equation obtained earlier in [17] in 
a different, diagrammatic way. However, instead of the zero- mode approx- 
imation made in [17], one can reduce this equation to the one-dimensional 
integral one, subject to a simple numerical analysis. 

This is done in section 3 of the paper. In case of two and more flavours 
one should use instead the bosonization method and find the chiral mass and 
chiral condensate from the steepest descent method, which will be a subject 
of the next paper. 

Finally in the last section a detailed discussion is given of the correspon- 
dence of our method with TLQ. 

2 Effective quark Lagrangians 

We consider in this letter the neutral dilute gas of N + instantons / and iV_ 
antiinstantons J, N + = AL = N/2, with the superposition principle [1-3] 



where 7 denote the set of collective coordinates for each instanton 7^ = 
{R^\ p^}- In what follows we for simplicity fix the instanton size = 
p and neglect perturbative gluon effects together with II, II interaction, 
assuming that the latter stabilizes gas at small density. 

Then the effective quark Lagrangian L E q L can be defined as 



N 



(i) 



i=i 





(2) 



where 




(3) 



i=i 




(4) 
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and 

N J4p(i) 

dl = IIW^) (5) 

i=l K 4 

The resulting L E q L is a polynomial of tp,^ of infinitely large power. 

To obtain L EQL one has to average in (2) over (5), 7 entering into A® 

as 

» [x,1) ~ Va ^ (x - my[{ x - my + P 2 ] w 

The straightforward use of the cluster expansion [18] for the averaging 
over yields 

Z = J D^Dij;e-^ f ^ s ^ dx+ ^i=^n =1 *«T> ( 7 ) 
where we have denoted (/, -flavour indices). 

= E flf = / da^J (8) 
/ fc =i 

and angular brackets imply integration over (Z7, while double brackets stand 
for cumulants [18] e.g. 

< 6» 2 >^R l =< ^ >C l R l - < >H l R I < ^ >^ lJ R l (9) 

Let us look at <C 9™ ^> more closely. The first term with n = 1 vanishes 
due to the integration over dfli. The next term with n = 2 will be of most 
interest to us in what follows. This can be written in the form (see Appendix 
1 for details of derivation) < 9f >= S/i/ 2 < ^/W 2 >? an d finally cast in the 
form, 

E < >= / ^r? r ^){(l " ^)X>^f (p)Jf(Pi) + 



- T^)E^°(py°(Pi) - 4E4? (p)Jf(Pi)- 

; E^(p)^i)} (10) 



2g 2 iV7 



4 



where q, d { are given in Appendix 1, and 



J?(p) = ^ + (p)t n O^(pi - q), Jr(pi) = V + (pi)i n O^(p + ?), 

jj»(p) = V> + (p)<rO^(Pi - 5), J?(pi) = ^ + (pi)# n O^(p + 5) (11) 

and n = F or 0, with F = 1, ...Nj — 1, and t° = l,t F is the flavour group 
generator. 

We have also defined in (10) 

(N?-1)V 4 ' 1 j 

and </?(g) is the Fourier transform of the instanton field (6) given in Appendix 
1, eq (A.6) . 

The 4q interaction term (10) is of the extended NJL type, containing all 
possible nonlocal bilinears. 

The advantage of (10) over the NJL Lagrangian is that it provides the 
in-built formfactor F(q) which makes loop corrections finite since it behaves 
as 

T<q) ~ V. (JV? -1)1^, , - < 13 » 

One can see in (13) that the small q behaviour of T(q) (which is due to 
the l/rr 3 asymptotics of the instanton vector potential(6)) implies massless 
exchanges between currents J(q) at large distances - a certainly unphysical 
feature, which is due to the ansatz (1) of noninteracting instanton gas. Taking 
into account the // and II interactions changes the instanton profile at large 
distances (possible forms are suggested in [3]) and makes T(q — > 0) finite in 
the realistic instanton gas or liquid model. 

There are several important properties of the new EQL (7). 

First of all, the sum over n is not limited, so therefore it contains vertices 
containing 4q,6q,... etc to infinity for any Nf. This differs from the standard 
ansatz of [9-12], where only the term with 2Nf quark operators appears. 

One can observe in (8) the contribution of all quark modes, since one can 
rewrite 9 l as follows 

6i = J dxi, + {x)A\x)^{x) = - E < > 4. < <l^o"V > (14) 

n.k 
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with e.g. < ip + \S 1 \u l n >= J ip + (x)(—id — im)u l n {x)d i x. and e l nk denned as 

e l nk = ~ < <|S A«Sok > (15) 

The standard TLQ prescription (9-12) obtains when one keeps in (14) 
only the zero mode coefficient, e l 00 , 

Consider now all terms in (7) with n > 2 and take the limit N c — > oo to 
perform the averaging over dfl explicitly. 

The EQL (7) is additive for instantons, therefore we consider it for the 
i-th instanton and omit the subscript i in what follows. One has 

OO 1 

n=2 U - i 

The double brackets in (16) imply subracting from < 9 n > disconnected 
pieces which appear while averaging in dQdR. 

Since each integration dR is accompanied by I/V4 the second term in the 
difference, 

< e n > R - < e k > R < e n - k > R , 

vanishes in the thermodynamical limit, iV — > 00, V4 
Hence the only terms left in 9 n ^> are of the form 

< e n > R ,n=« n >n-< e 2 > n < B n ~ 2 > Q 

Now in averaging over dVL one obtains in the leading terms of the 1/N C expan- 
sion the set of disconnected (in color indices) products, which are subtracted 
in (18) and (n — 1)! equivalent connected terms, which differ only by permu- 
tation of indices. Thus one comes to the EQL defined in (16) 

x (27r)V(E^)^(n n^A-i + i 1 - 1) (19) 

i=i fc=i i=i 

The EQL (7) and (19) is the central point of this letter and in the following 
section we use it to derive equations for the chiral mass. 



(17) 

— /•/ 

v 4 



00, 7T = const. 



.... >R (18) 
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3 Chiral mass, chiral symmetry breaking and 
the gap equation 

The multiquark terms <C 6 n ^> in the effective Lagrangian (7) contain all 
information about possible quark mass creation which automatically signifies 
CSB. The form of < 9 2 > in (10) tells us that the averaged 4q interaction is 
of the extended NJL type [15] and is nonlocal. 

Nevertheless one can use for the 4q term the same gap equation, derived 
from the Dyson-Schwinger method, as it was done in the NJL formalism 
[15]. However the essential difference of our expression (7) from the NJL 
Lagrangian is that the former contains all powers n, i.e. all n-quark effective 
vertices, n — 4, 6, ...oo . 

Therefore the gap equation, defining the creation of the chiral mass in 
our case, should contain in general this infinite number of vertices, unless 
one can find a small parameter suppressing terms with large n. To answer 
this question we shall consider first the case of one flavour and compute the 
averaged quark propagator in the instanton gas. 

We shall derive it in the large N c limit and shall obtain the integral 
equation for the chiral mass, containing contributions of all n-quark vertices. 

At this point one can find from (19) the effective quark mass in the 
mean-field (Hartree) approximation. To this end one should identify in L 
(19) the structure —% J j0^'4 ,+ {p)M{p)il){p). In the leading order at large 
N c one replaces in (19) (n — 1) colorless bifermion combinations by their 
vacuum expectation values, and in Hartree approximation by the effective 
quark propagator S(p) 

S(p) = - \- ( . (20) 
p-iM(p) 

Thus each of the colorless bifermions yields a factor N c and we obtain 

at oo „ n j n 
4 c n=2 k=l v / i=l 

+ (21) 

where the last bracket in (21) refers to the contribution of antiinstantons. 
One can rewrite this expression in a simpler form 
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M(p) = +ij±y £ < P\(S ~ (ATT'IP >R (22) 

Eq.(22) is equivalent to Eq.(7) of [17] when one does Fourier transforma- 
tion and averaging over dRi. 

Since S(p) contains M(p) as in (20), two coupled equations (20) and (21) 
define M(p) through A(p) and the density parameter jp^- 

One can introduce the "Green's functions" G(G) for the quark propaga- 
tion in the in the instanton (antiinstanton) field A(A) with effective mass 
M(p), e.g. 

(P _ lM _ A)G(p, p') = (2tt) 4 5( P - p>) (23) 
Separating 5-function term from G, one obtains 

(2^-^) ^/) m 

^ p >~ p-iM{p) ^ (p-iM(p))(p> -iM(p>)) { ' 

and similarly for G and £. Insertion of (24) into (22) yields 

+ (A^A,(^£)} (25) 
and for £ one has an equation following from (24) 

= Mp - P ') + /Mp- q)j -± m((q , P >) (26) 

and similar equation for £(p,p') with replaced by the antiinstanton field 

Since one expects M(0) to be around 200-300 MeV, [9-12, 17], the spectral 
expansion of G and the ZMA done in [17] which is valid for vanishing M(0) 
is not convincing. Therefore one has to solve Eqs. (25-26) numerically. 

The calculations done above in this section, demonstrate that (19) is the 
EQL corresponding to the Pobylitsa equation [17] obtained for Nf — 1 and 
N c — > oo. In case of Nf > 1 the EQL is given in (7) and its explicit form for 
n = 2 in (10). For many flavours (Nf > 2) it is natural and (easier) to use 
the bosonized form of EQL, which will be published elsewhere. 
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Finally we discuss in the section symmetries of our EQL, J2n < , > . We 
consider only the case N + = AL, so that the net topological charge of the 
system is zero, and there are no global zero modes. Correspondingly the 
effective 'tHooft Lagrangian [7] with the global zero modes does not occur 
in this case. Now the 4q term in (7) can be rewritten as in (A. 9), the form 
which respects U(Nf) x U(Nf) symmetry. The same can be shown for all 
terms in (7), and therefore violation of SU(Nf) x SU(Nf) can occur only 
spontaneously while C/a(1) has anomaly. The explicit example of CSB is 
given by a nonzero solution of (25-26), derived above in the limit N c — > oo. 

4 Summary and discussion 

The main result of the letter is the EQL (7), with explicit form for n = 2 
and any Nf in (10), and for all n and Nf — 1, N c — > oo in (19). 

The 'tHooft-type ansatz of EQL of ref. [10-12] is obtained from our EQL 
(7) when one selects there n = Nf and keeps in the mode expansion (14) only 
the zero mode contribution. One can check, that in this case automatically 
appears the determinantal structure of the 'tHooft ansatz. We have not 
found any arguments why nonzero modes in (14) and all terms with n ^ Nf 
should be dropped. 

It is interesting, that the structure of the 4q term in the EQL Eq. (10), is 
so much resembling the extended NJL model with all Lorenz group channels 
(see (A. 12)), in contrast to the standard TLQ Lagrangian [10-12], where in 
the leading 1/N C order only scalar and pseudoscalar channels appear. 

It is not impossible, however, that the resulting effect of all terms of EQL 
in the sum over n in (7), will be numerically close to the standard TQL results 
[10-12] for M(p) and < qq >. As one example of this sort one may quote 
the equation (22) for the effective quark mass which is obtained above from 
our EQL in the Hartree approximation and which coincides with that of ref. 
[17], where diagrammatic analysis for the quark propagator was exploited. 

The author of [17] applied the ZMA for (22), valid in the limit of small 
M(p), and obtained a solution for M(p) closely resembling that of the stan- 
dard TQL analysis [9-12]. However the resulting value of M(0) ~ 350MeV is 
not small and this fact casts some doubts on the whole numerical procedure 
of [17]. The crucial check of this coincidence would be the exact solution of 
Eqs. (25-26) for M(p) and its comparison with the standard TQL solution of 
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M(p) [9-12]. This will be a subject of a future publication. 

As it is, our EQL (7) is an exact result for any N c and any Nf, which 
enables one to calculate different physical effects, from qq creation due to 
the 4q vertex, to chiral symmetry breaking effects (M(p), < qq > etc). The 
corresponding detailed calculations and comparison with the existing TQL 
results is planned for the future and will help to clarify the mechanism of 
chiral dynamics in the QCD vacuum. 

The author is grateful to G.'tHooft for a valuable discussion, to D.I.Diakonov 
for a useful discussion and correspondence, and to B.O.Kerbikov and D.S.Kuzmenko 
for numerous discussions in the course of the present work. 

This work was supported in part by grants RFFI 96-02-19184, RFFI-DFG 
96-02-000 88G and INTAS 94-2851. 
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APPENDIX 1 



Derivation of the 4q term of the effective quark Lagrangian 



< 6 2 >=< J ^(x)A(x)^(x) J ^(y)A(y)^(y)dxdy > Q , R = 

J (^ h (xh^(x))dx J (^%yh,M 2 (y)dy \^f (A.l) 
where a,/3 are color indices, and 

2 r 

x a(x - R)a{y - R) (A.2) 

where 

a(a;) = ^T7) (A ' 3) 

and the sign +(— ) refer to the instanton (antiinstanton) case respectively. 

In (A.l) one can notice that color indices are intercharged in both quark 
bilinears, which calls for the application of Fierz transformation, which gives 

1 1 

7m x 7 „ = 1 x 1 - -( 7/U x 7 „) - -( 757 ^ x 75 7 M ) - 7s x 7s (A.4) 
Note that we are using Euclidean 7 matrices, i.e. 

74 = ( 1 _ x ), 7* = 75 = ( _j 1 ), 7^ + 7^7m = 2^ 

For the following it is useful to go over to the momentum space and to 
this end we need the Fourier transform of the instanton vector potential 



where 



A, a (x -R) = 2% auP W J J0p . qMq)f M*-m (A.5) 
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< 02 >= S ^^^^^ h ^^m h ^h^ 



Insertion of (A. 5) into (A.l) yields 

rV) 2 

(27r) 4 (27r) 4 27r) 4 (A r c 2 -l)^4' 

(^ii^Q 2 - QuiQuJ^io) ( A - 7 ) 
Using (A. 4) and a similar Fierz relation for scalars 

1 1 

1 x 1 = -{1 x 1 + 7m x 7m + -a^ x - 75 7 M x 757^ + 75 x 75} (A.8) 

one obtains for (A. 7) 

^2„2\2,„2 



< 



1=1 

(A.9) 

where Oi and q, are 

Oi = l,i^ 7s7a» 7s, <V = liPlv 2 - lvlii for i = 1,2,3,4,5 

Q = 1,-^,-^,-1,0 (A.10) 

, _ 1 1 111 

c * ~ 4'4'~4'4' 8 

Now we take care of the flavour indices 

E = -o4r^^ + (A - n) 

to rewrite the curly brackets in (A.9) as 

{...} = 2q 2 £ c^ + (p)0^(Pi - qm + (pi)0, t t F ^(p + q))+ 

i,F 

+JT E c t ^ + (p)O t l F ^(pi - q))^ + ( Pl )O t l F ^(p + q))+ 

f i 



12 



-2 £ cj^ + (p)g0^(pi -?))ty + (Pi)50i* J V(p + ?))- (A.12) 

F,i=l 

^/ i=i 

We now can introduce currents 



^ ,iy (p) = ^ + (p)Oi(^,i*)^(Pi-?) 

r for j/ 

in (A.12) as 



and similarly for J i ' , where instead of Oj stands gOj, to rewrite all terms 



{...} = 2g 2 $>jf (p)jf + ^£qJ°(p)J°( Pl ) 



-2E4Jf(p)Jf(pi)-^E^^(p)^(Pi) (A.13) 

i,F ly f i 

As a result (A. 9) assumes the form 

2 y ^ Pl dg2(47rV)V 2 (g) 2 1 F F 



< 



#(1 - 4)E^°(p)^(Pi) - 2 £^f (p)Jf (pi)- 

JV / JV c i i,F 

-lJ-E^(P)-?(Pi)}- ( A - 14 ) 
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